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Abstract. We examine bootstrap percolation on a regular {b + l)-ary tree with initial 
law given by Bernoulli(p). The sites are updated according to the usual rule: a vacant site 
becomes occupied if it has at least 9 occupied neighbors, occupied sites remain occupied 
forever. It is known that, when b > 6 > 2, the limiting density q = q{p) of occupied 
sites exhibits a jump at some pj = pj{b,9) 6 (0,1) from qj := qipj) < 1 to q{p) = 1 
when p > Pj. We investigate the metastable behavior associated with this transition. 
Explicitly, we pick p = pj + h with h > and show that, as h I 0, the system lingers 
around the "critical" state for time order /i^^^^ and then passes to fully occupied state in 
time 0(1). The law of the entire configuration observed when the occupation density is 
q 6 {qj, 1) converges, as h J, 0, to a well-defined measure. 



1. Introduction 

Consider a connected and bounded-degree graph G = {V,E). Bootstrap percolation 
on G is a stochastic particle system on the configuration space {0, 1}^ defined as follows: 
Given a configuration cr = (ax) G {0,1}^, we interpret = 1 as "x is occupied" and 
cr^ = as "x is vacant." Pick an integer > 1 to be referred to as the threshold. Under 
continuous dynamics, a vacant vertex becomes occupied at rate 1 as soon as it has at least 
6 occupied neighbors, otherwise the rate is zero. An occupied vertex remains occupied 
forever. In discrete-time evolution the updates are performed at all vertices simultane- 
ously at regular time intervals (which means that the evolution is entirely deterministic 
and the model is thus an example of a cellular automaton). The initial configuration is 
sampled from an i.i.d. distribution; i.e., Bernoulli percolation with parameter p. Hence 
also the name of the process. 

Bootstrap percolation arose in physics context in studying the effects of dilution on 
magnetic properties of crystalline materials (Chalupa, Leath and Reich 0). Notwith- 
standing, the main source of motivation for probabilists was the connection with metasta- 
bility, and generally non-equilibrium statistical mechanics (Aizenman and Lebowitz [IJ). 
A large body of work exists for the problem on the hypercubic lattice G = Z**. In this 
case, the dynamics converges to fully occupied lattice whenever p > and 9 < d, and 
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never to fully occupied lattice when p < 1 and 6 > d (Schonmann fr2\ ). While the "crit- 
ical density" in the most basic case d = 2 = 9 is zero (van Enter IHl), the transition from 
partial to full occupancy is exhibited in the limit p J, by the process restricted to boxes 
of scale exp('^/p) as A varies through ^^/ig (see Aizenman and Lebowitz [IJ, Holroyd IITOl 
for results in d = 2 and those by Cerf and Cirillo [5J, Cerf and Manzo IH, Holroyd llTT| 
and Balogh, Bollobas and Morris 121 dealing with higher dimensions). 

In this short note we examine bootstrap percolation for G being a homogeneous tree 
with degree b + 1. The dynamics, both discrete and continuous-time, can in this case 
be characterized rather explicitly. (This is not always the case even on a tree; e.g., the 
stochastic Ising model or kinetically constrained models of Ising type.) It is well known 
(Chalupa, Leath and Reich [7], Balogh, Peres and Pete [3J) that, for & > > 2, there 
exists a critical density px such that the limiting density of occupied sites converges to 
c]j = cjj{p) < 1 when p < pi and to one when p > pj. There is yet another critical value, 
Pc with pc < Px, above which the terminal configurations will contain infinite connected 
components (Pontes and Schonmann [9J). 

Our main goal is to point out an interesting metastability effect associated with the 
first-order transition at p = px. Our findings are as follows: Por the initial density p 
slightly above px, i.e., p = pj + h, there is a metastable time scale of length 

t{h):=och-^/^ + 0{l), hiO, (1.1) 

for some a G (0, oo), at which we see a transition from density qj to full occupancy in a 
window of times of order unity! This extremely sharp "cutoff phenomenon" is perhaps 
even less intuitive when we add that, as h [ 0, the law of the full configuration at the 
time when density equals some value q S {qj,'^) tends to a distinct measure }iq (of that 
density). Thus, a few moments before t{h) we see a distinctly non-fully occupied graph 
while a few moments after t{h) the graph is essentially full. Presumably, the sharpness 
of the flip from the metastable to the stable equilibrium is a result of some global event 
happening at spatial distance — measured in graph-theoretical sense — of order h"^^'^. 
However, the nature of this event remains mysterious at this point and is definitely 
worthy of further investigation. 

We finish by noting that analogous first order transitions to the one studied here have 
been found of relevance for jamming (Schwarz, Liu and Chayes [13]). A somewhat 
related model mimicking a dynamic of single avalanche in directed tree geometry has 
been thoroughly analyzed by Biskup, Blanchard, Chayes, Gandolfo and Kriiger ||4|. 

2. Recursion equations 

We begin by presenting the basic recursion-equation calculations that are permitted by 
the tree structure of the problem. We will examine two distinct dynamics; first, a discrete 
time process, when at integer times each unoccupied vertex examines its neighbors, and 
if there are more than Q of those occupied, it becomes occupied. The second, continu- 
ous time, is such that each vertex is equipped with a clock that rings at exponentially 
distributed (independent) time intervals, and the update happens accordingly. 
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We will consider the process on the {b + l)-ary regular tree T^, and denote the con- 
figuration at time f by (7"(i) = (c^j(O)xeTj- Here ax{t) G {0, 1}, with indicating vacant 
and 1 occupied state, and the initial condition a{0) is i.i.d. with density p. The char- 
acter of the dynamics implies that f i— > crx{t) is increasing at each x, and can thus be 
completely characterized by the quantities 

T:, ■.= mf{t>0: a^it) = 1}. (2.1) 

In particular, T;,; = if x is occupied initially. Define 

Pp{t) := P{T, < t) (2.2) 

which, by the symmetries of the dynamics and initial data with respect to tree automor- 
phisms, does not depend on x. Let Bin(b, q, 9) denote the probability that a sum of b i.i.d. 
Bernoulli((j) is at least 9, i.e., 

Bm{hq,9):=p^(J^q\l-qf-K (2.3) 

Our first observation is: 

Lemma 2.1 Let Z = Ifor the discrete-time process and Z = exp(l) /or the continuous-time 
process (independent of the a's). Then 

Pp(0 = p + (l-p)£(Bin(& + l,Qp(t-Z),0)l|z<f}), (2.4) 
where t ^ Qp{t) is the unique solution in [0, 1] of 

Qp(0 = p + (l-p)£(Bin(&,Qp(i-Z),0)l{z<t})- (2.5) 
In particular, for the continuous time evolution 

AQp(i) = Wp(Qp(0) (2-6) 

where 

Wp{q) := p + (1 - p)Bin{b,q,9) - q. (2.7) 

Proof. Let denote the rooted tree with forward branching number b and let denote 
the root. Consider the bootstrap process restricted to T^. We claim that the quantity 

Qp(0 := P(T0 < t) (2.8) 

satisfies ( |2.5| |. Indeed, let S denote the first time the root has at least 9 occupied neighbors 
and let t — Z be the time of the most recent "clock-ring" at the root before time t. Then 

{0 < T0 < = {T0 > 0} n {S < T0} n {S < i - Z}. (2.9) 

But on the event {S < T^} each of the neighbors at the origin evolves independently 
according to same law as the origin itself. Hence, to get the event on the right-hand side, 
we need that at least 9 neighbors of the origin are occupied by time t — Z, i.e., 

P(0< T0 <f) = (l-p)£(Bin(&,Qp(i-Z),0)l|z<f}), (2.10) 
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where the prefactor accounts for the probabiHty that, at time zero, the origin was vacant. 
This immediately yields ( |2.5[ ). 

Now let us go back to and the full bootstrap percolation process on this graph. 
Let X S Tf, be any vertex. If Sx is the first time x has at least 6 occupied neighbors, 
then {0 < Tx < t} can again be written as > 0} intersected by {S^ < T0} and 
{Sx < t — Z}. The difference compared to the "directed problem" is that now x has b + 1 
neighbors and so 

P(0 <Tx<t) = {l- p)£(Bin(& + l,Qp{t - Z),Q) l|z<t}). (2.11) 



This again immediately gives ( |2.4| |. 

To get ( |2.6| l we note that, in the case of continuous-time dynamics, Z is exponentially 
distributed. A simple change of variables makes ( |2.5| l into 



Qp(0 = p+(l-p)^'dze-('-^)Bin(&,Qp(z),0) (2.12) 

This shows that i 1— > Qp(i) is smooth; differentiation then yields ( |2.6| |. □ 
Remark 2.2 For discrete-time evolution, the analogue of ( |2.6| l is 

Qp{n + l)-Qp{n) = Wp{Qp{n)) (2.13) 

and so, by standard approximation methods, the discrete and continuous problems can 
be shown to follow similar trajectories. However, to avoid tedious calculations, we will 
not attempt to cast this statement in quantitative form. 

Since Qp{t) readily determines Pp{t), we will focus the forthcoming discussion on Qp. 
The indisputable advantage of continuous time is that ( |2.6| can be integrated out: 

(214) 

The analysis of Qp{t) thus reduces to the analysis of the denominator Wp. Here we 
notice the following facts (some of which can be traced to calculations in Pontes and 
Schonmann ||9| and earlier work on this model): 

Lemma 2.3 Suppose b > 9 >2 and < p < 1. 

(1) There is a q G (0,1), independent of p, such that q 1-^ ^p{q) is convex for q < q and 
concave for q > q. 

(2) There is p G (0, 1) such that q ^^piq) is decreasing on [0, l]for p > p (with a unique 
local minimum at q = 1) while for p < p, there is a qmm{p) S (0, 1) such that the local 
minima of q ^ Wp{q) in [0, 1] occur at q ^ {qminip), !}• Both of these are then strict. 

(3) p ^ Wp((^jnin(p)) is increasing, so there is pj G (0, p) such that Wp{qmm{p)) < ^ for 
p <Pt and Wp(^mm(p)) >Oforp > pj. 

(4) When Wp{qmm{p)) < 0, then q^mip) > P- 
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Proof. To get (1) we notice that, by a coupling of Bernoulli's to independent Uniform([0,l]) 
random variables, 

^Bin(&, q,e) = b(^^l (1 - qf-'. (2.15) 

For b > 9 > 1, this starts increasing and then becomes decreasing, implying the said 
convexity types of Wp. Note that Wp(0),Wp(l) < and set a := ^Bm{b,q,0). If 
Wp{q) = fl(l — p) — 1 < then q Wp{q) is (strictly) decreasing throughout [0, 1], while 
for fl(l — p) > 1, there is a secondary (strict) local minimum in (0, q). This yields (2) with 
p defined by fl(l — p) = 1. 

To get (3), we note that p ^ Wp{qram{p)) has derivative 1 — Bm{b,qram{p),0) which 
is positive for all < p < 1. The crossing of zero level must occur as, on the one hand, 
Wp(0) = p while Wp(0) = —1 and so q ^ ^p{q) definitely plunges below zero for 
small p, while, on the other hand, Wp is strictly positive on [0,1) for p close to p. To 
see (4), we note that the (strict) convexity on {0,q) implies 

Wp{q)>Wp{0)+W'p{0)q = p-q (2.16) 
and so Wp{q) > ior q < p. □ 

Remark 2.4 We note that q ^p{q) is convex throughout [0, 1] when 6 = 1, while it is 
concave when 6 = b. No local minimum thus appears inside (0, 1) for either of these 
cases and Wp is thus minimized (and vanishes) only at q = 1. None of the interesting 
effects described here occur in either of these cases. 



The findings in Lemma 2.3 imply that that, for p < px — which by (4) forces p < 
qminip) — as i ^ 00, Qp{t) tends to the smallest positive root of Wp while for p > pj, the 
integration eventually passes through the bottleneck at qmm{p) under the graph of Wp 
and so Qp{t) ^ 1. A more detailed look at the bottleneck then yields: 

Theorem 2.5 Suppose b > 6 > 2. Let px denote the value of p at which min Wp = 0. Set 
qT '■= qminipT) dud let p;, := px + h. Denote 



2"pt 

Then for each r G IR, the limit 



, " J\^. (2.17) 



cpir):=\\mQp,X^h-'^' + r) (2.18) 

exists and is a solution to the ODE 

cp'{r) = Wp^{cp{r)), rGR, (2.19) 

subject to boundary conditions 

lim (p{r) = qj and lim (j){r) = 1. (2.20) 
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(This solution is determined uniquely up to a shift of the argument.) Equivalently, given q G 
(qjA), ifhiq) is defined by Qp^{th{q)) = q, then 

thiq) - och-^'^ (p'\q), qT<q<l. (2.21) 

hlO 

This result encompasses the principal findings of this paper. The interpretation is 
that, for initial density p slightly above pj, the system evolves into a metastable state 
close to that for the critical initial density pj. However, this metastable regime ends 
rather abruptly at time K{p — pj)^^^'^ + 0{1) when a sudden transition to full occupancy 
occurs in a time window of order unity. 



3. Proof of metastability phenomenon 



lim 

010 



We now commence with the proof of Theorem 2.5 Key in the calculation is the following 
observation: 

Lemma 3.1 Let q < 1 and let wbea non-negative, fiinction such that 

w{q)=0 = w'{q) and w"{q)>0. (3.1) 
Then for all sufficiently small 3 > 0, 

f^^^ dx 1 

/ — ^7 r — 59^^^'^ exists and is finite, (3.2) 

Jq-s iv{x) + 9{l-x) ^ \ ' ^ ' 

where j6 := n/ ^J^w"{q){l — q). 

Proof. This follows by standard, albeit tedious, calculations. First, set x q + {1 — q)x 
and w{x) = (1 — q)w{x) and note that this reduces the problem to q = 0. (We also need 
to note that w"{q) {1 — q) = w"{0) and that 3 gets to be scaled by the factor (1 — q).) 

We will henceforth assume that q = and that S is so small that w{x) > for all x G 
[— ^, \ {0}. First, abbreviate 

1 1 

^^^^ w{x)+e{l-x) ~ w{x) + 9 ^^'^^ 

and check that 

R{x) + R{-x) 

^ . [w{-x)-w{x)][w{x) + w{-x) +19^x9] 

[w{x) + 9] [w{x) + 0(1 - x)] [w{-x) + e] [w{-x) + 0(1 + x)] ■ ^ ' 

The bounds w{—x) — w{x) = 0{x^) and x^/w{x) = 0(1) show that this expression is 
bounded uniformly as J, and |x| small; the J, limit is then obviously zero. In 
light of the X — X symmetry of the integration domain and the Bounded Convergence 
Theorem, we may thus replace the integrand in ( |3.2[ ) by [w{x) + 9]^^; the claim thus boils 
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down to showing that 



lim 



dx 



Here we abbreviate 



and invoke the identity 



s w{x) + 9 
1 



fie 



-1/2 



exists and is finite. 



1 



fl2 := -w"{0) and fls := -w'"{0) 



(3.5) 



(3.6) 



1 



+ 



2w(x) + e 2w{-x)+e a2X^ + e 



w{—x) — w{x) 



[a2X^ + e] [w{x) + e] [w{-x) + 6] 



1 w{x) — UzX^ — a^X^ 1 w{—x) — fl2^^ + ^3^'^ 



2[a;(x) + 0][fl2x2 + 0] 2[w{-x) + e][a2X^ + 0] ^^'^^ 

Again, all three expressions on the right-hand side are bounded uniformly in x G [— (J, ^5] 
and have a well-defined limit as J, 0. The Bounded Convergence Theorem again en- 
sures that the integral of [w{x) + 6]^^ equals that of [a2X^ + 9]^^ plus a quantity that has 
a finite limit as J, 0. The claim now follows by noting that 

dx Tz 2 



-s fl2^^ + 9 ^a29 
which is checked by a direct calculation. 
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+ 1, 



no. 



(3.8) 

□ 



Proof of Theorem^Sj It suffices to prove ( |2.21[ ); indeed, the claim ( |2.18 1 then follows by the 
strict monotonicity of t Qp(0 ^rid of the solution to the ODE ( 2.19| subject to ( |2.20[ ). 
First note that 

ds 

hici) = I ^X7^. (3.9) 



Thus, the point- wise limit Wp,^(s) Wp^(s) and the fact that Wp^(s) > for s > qj 
imply that t]j{q) — t]t{q ' ), for q,q' > qj, converges to (p^^{q) — (p^^{q'), for some solu- 
tion r I— > (p{r) to ( 2.19| -2.20 1. To show ( 2.21| | it now suffices to prove that (for at least 
one q elR) 

lim\ th(q) — ah^^^^] exists and is finite. (3.10) 



A calculation yields 



Ph 



;wp^(s) + 0,(i-s)) 



(3.11) 



1 - Pt 

with 9}j := h / [1 — p]j) . Noting that Wp^ is uniformly positive on the complement of the 
interval {qj — S,qT + S) and observing the asymptotics 

l^-l = o{h'/') and ^ 
l-px h 



(1 - px) = o{h'% 



(3.12) 
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to get ( 3.10| l it thus suffices to show that 

r fn+^ ds -1/2I 
lim / — — ^-T — ^ — B6 ' exists and is finite (3.13) 

eioU^j-s Wp^(s) + 0(1 -s) ^ J ^ ' 

for j6 := n/ ^JlW'^^{qT){l - qj). But this follows from Lemma [sjj □ 

Remark 3.2 The asymptotic analysis in the above proofs allows us to control the speed 
of convergence at the times way before (or after) ah^^^^ -\- 0{1). Thus one can prove, for 
instance, that for all < A < a, 

h-'/^[Qp,X^h-'^^)-qj] ^Citan(c2A-C3) (3.14) 
for some constants Ci, C2 > and C3 G R. 



4. Process convergence 

The above calculations give us information about the evolution of the law at a single 
vertex. However, as we will show in this section, the tree structure and some additional 
arguments allow us to control the law of the entire infinite-volume configuration at times 
close to and after passing through the metastable regime. For f > 0, let jip^t to denote the 
law of {ax{t)) started from Bemoulli(p). 

Theorem 4.1 There is a stochastically increasing family of tree-autormorphism invariant mea- 
sures {vq)o<q<i on {0,1}^* with density Vq{crx = 1) = Bin(£' + l,q,6) such that: 

(1) iff — Pt dnd q := limf^oo Qp{i)r then 

^P,t -^Vq, t ^ 00. (4.1) 

(2) lfp>pi then 

}ip,t vi ■■= h, t 00. (4.2) 

In addition, we also have: 

(3) Set Pij = p-Y + h and, given q G {qj, 1), let th{q) be defined by Qp,^ {thi'])) = ^- Then 

Remark 4.2 We note that the statement (3) is valid only for the continuous-time dynam- 
ics as, for discrete time, even the occupation density changes in discrete quanta. 

The convergence in (1-2) follows trivially from the fact that t 1— > jip^t is stochastically 
increasing; the stochastic ordering oi q ^ Vq is a consequence of the existence of mono- 
tone coupling of the processes with different initial densities and — in part (3) — also 
the monotonicity of the dynamics and of the function q 1— > tfi{q). The main issue is thus 
the proof of (3). The key observation is the fact that the effect of the rest of the system on 
a finite set can be represented through a time-dependent boundary condition. 

Let A be a finite connected subset of and let dA denote the vertices in Tf, \ A that 
have a neighbor in A. For each x G dA let T(x) denote the (unique) subtree of Tj, \ A 
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rooted in x. For any initial condition (Tx{Q) = (t^, let {(Tx{t): x G T;,} be the state of 
bootstrap percolation on Tj, at time t. In addition, for each x G dA, let {yjA,z{i) '■ z G 
T(x)} be the state at time t of bootstrap percolation on the graph T(x) which is started 
from the same initial condition as a. Given a trajectory {rjA,z{i)'- z G dA, t > 0}, we 
now define {f/A,.T(0 : ^ G ^} to be the state of bootstrap percolation on A at time t that 
initiates from rjA,x = o'x and evolves against the time-dependent "boundary condition" 
Va,z{^) at allz G dA. 

Lemma 4.3 The laws of {crx{t) : x ^ A, t > 0} and {r]A,x{i) x ^ A, t > 0} are identical. 

Proof. For each x G dA, let n{x) G A be the unique neighbor of x in A. We can now 
couple the evolutions of {cxit) : x G Tj,, t > 0} and '■ ^ £ T^,, f > 0} by using 

the same rings of Poisson clocks. Given a sample of all these rings on Ty, for all x G dA 
we have 

(^z{t) = r]A,z{t), z G T(x), < f < T„(,,). (4.4) 

However, at i = T^j-;,,) the vertex n{x) becomes occupied and, since T(x) has no other 
neighbor in A, the evolution on T(x) becomes irrelevant for that in A. It follows that 
C7z(f) = f]A,z{i) for all z G A and all f > which implies the claim. □ 



As we will see. Theorem 4.1 will be a consequence of this, slightly stronger, claim: 



Proposition 4.4 Let a. he as in { 2.17) and let h ^ r^he such that both rf, and ah _ j,^ 



tend to infinity ash [Q. Define G {0, 1} x R via 

Sx:=(l{T.<o,}^7^.-«^"'''l{r.>,,,})- (4-5) 

Ash [Q, the law of [Sx) induced by Pp with p = pj + h on the product space ({0, 1} x M)'"^* 
tends to a probability measure which is independent of the choice ofrf^. 

Proof. We begin by noting that the probability of the event 

{Tx <ri,}U {ah-^'^ -rh<Tx< oih'^'^ + (4.6) 

tends to one as h J. 0, for any r;, as above. Indeed, the events are disjoint once h is 
sufficiently small, the probability of the former is at least q-x in the h [0 limit by the fact 



that Qpj+h{i) > Qpt(0 ~^ '?T while the probability of the latter tends to 1 — by ( |2.21[ ). 
It follows that the cutoff to Tx > r^^ may be replaced by 

_ < 7^ < + (4 7) 

in the second half of ( |4.5| . 

To describe the convergence of the full distribution, pick a finite set A C Tj, and, 
given the initial data a a '■= {cTx'- ^ £ A} and boundary data S^a '■= {Sx - x G 9A}, let 
Ha{ — \ SdA, cta ) be the bootstrap percolation process on A started from the initial configu- 



ration a A and evolving against the time-dependent boundary data S^a- Lemma 4.3 tells 
us that the law of the variables {Sx'. x G A} is that of }iA{ — \5dA,o'A) integrated over 
Bernoulli(p) variables aA and i.i.d. boundary data where each Sx, x G dA, is sampled 
independently from the law of at the root of a rooted tree. 
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Recall that, by { 2.21) the law of S^/^ has a non-degenerate weak limit as ?2 J, 0, which is 
the same regardless of how r;, oo and ah^^^'^ — ^ co. By definition, also the law of 
aA trivially converges ash [ 0. As there is only a finite number of possible (Ta's, to prove 
convergence it suffices to show that, for each a^, the family of measure- valued functions 

SdA^ HA{-\SdA,0-A) (4.8) 

is equicontinuous. For this we note that a change in Tx for x G dA by e > will only 
be felt inside A if the Poisson clock at the (unique) neighbor y G A of x clicks within 
the interval [Tx, Tx + e]. This has probability of order e. The functions ( |4.8| are thus 
Lipschitz continuous in the boundary values and convergence follows. □ 



Proof of Theorem 4.1 '3). Assuming r^ < ti,, we have 



l{T-x<i-/,} + ^{T,>rn} l{T^<f,,} • 



The convergence of f/p,„t;, then follows from Proposition 4.4 and ( 2.21\ 



(4.9) 

□ 
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